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Direction-dependent anisotropic exchange is a common feature of magnetic systems with strong 
spin-orbit coupling. Here we study the effect of such exchange upon macroscopic magnetic anisotropy 
for a face-centered-cubic model. By several theoretical techniques, we show that, both in the para¬ 
magnetic and ordered phases, the magnetic anisotropy is induced by fluctuations. Moreover, the 
magnetic anisotropy differs in the paramagnetic and ordered phases: in the paramagnetic phase 
the susceptibility is maximum along the (111) directions, while the magnetic moments orient along 
(110) or (100) in the ordered phase. We suggest that such “anisotropy switching” can be a common 
feature of strongly spin-orbit coupled magnets. 


PACS numbers: 75.10.Jm, 75.30.Gw, 75.10.-b 

Insulating magnetic compounds with heavy elements, 
such as Ad and 5d transition metal ions, are of consid¬ 
erable interest in the study of quantum magnetism. In 
these materials, strong electron correlation and spin-orbit 
coupling often lift the orbital degrees of freedom, realiz¬ 
ing “pure” spin systems with effective pseudo-spins com¬ 
prised of by entangled microscopic spins and orbitals. 
Generically, such strong spin-orbit coupling gives rise 
to direction-dependent exchange interactions. The pres¬ 
ence of such interactions potentially affects the magnetic 
properties in a significant way, possibly stabilizing ex¬ 
otic phases, such as spin-liquid statesExperimentally, 
several iridates^'— and double perovskitesi^riS have been 
studied in these contexts. 

In this paper, we theoretically study a model of an an¬ 
tiferromagnet on a face-centered cubic (fee) lattice with 
direction-dependent interactions. This model may apply 
directly to cubic double perovskites, but also serves as a 
representative example of the general class of frustrated 
systems with directional couplings. We observe that, at 
low temperature, the primary role of the directional in¬ 
teractions is to partially lift the ground state degeneracy 
of the isotropic problem, selecting a particular subset of 
ground states. In addition to this “direct” lifting of the 
ground state manifold, quantum fluctuations also affect 
the ground state selection - a phenomena known as “or¬ 
der by disorder” 

At the same time, the directional couplings and fluctu¬ 
ations also affect the magnetism in the high temperature 
paramagnetic phase, and in particular the magnetic sus¬ 
ceptibility and anisotropy. Importantly, the nature of 
fluctuations in the high temperature phase is different 
from those in the ordered phase at low temperature. In 
the high temperature phase, since the spins are strongly 
disordered, correlations between nearby spins are domi¬ 
nant. On the contrary, long wavelength fluctuations tend 
to be more important in ordered phases. Hence, the con¬ 
tribution of fluctuations in these two regimes may result 
in qualitatively different behavior. 

We show that this is indeed the case in fee antifer- 


romagnets. By a combination of the Luttinger-Tisza 
method, spin wave theory, and the high-temperature 
expansion, we show that the magnetic anisotropy 
“switches” between the paramagnetic and the ordered 
phases. Specifically, at low temperature in the ordered 
phase, either the ( 100 ) or ( 110 ) directions are favored, de¬ 
pending on the sign of the directional exchange coupling, 
while in the paramagnetic phase, the ( 111 ) direction is 
favored. We show that quantum fluctuations dominate 
in selecting the ground state, while thermal fluctuations 
drive the magnetic anisotropy in the paramagnetic phase. 

The model we consider consists of spins on the fee lat¬ 
tice with direction-dependent nearest-neighbor (NN) ex¬ 
change and isotropic second-neighbor Heisenberg inter¬ 
action, 

H = Ho+Hh ( 1 ) 

with 

Fo = -J2^Si-S,- 

ihj) [hj] 

( 2 ) 

ihj) 

= (3) 

i 

Here, = (5'f, Sy, Sf) is the spin operator for the local¬ 
ized spin at ith site and Sij is the vector connecting ith 
and jth sites. We choose the standard convention so that 
the conventional unit cell of the fee cube has unit length. 
The first and the third sum in Eq. m is taken over all 
the NN sites and the second one is for second-neighbor 
sites. Ji and J 2 are the NN and second-neighbor Heisen¬ 
berg interaction, and J 3 is the direction dependent in¬ 
teraction. Hy is the coupling of uniform magnetic field 
h = {hx, hy,hz) to the spins. In the following, we focus 
on the case Ji, J 2 > 0 . 

Classical ground states: We first consider the magnetic 
ground states at h = 0. When J 3 = 0, the classical 
ground state of the model in Eq. (HD is given by a four 
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FIG. 1. Schematic picture of the classical ground state for the model in Eq. [T]with (a) J 3 > 0, and (b) double-q and (c) 
triple-q magnetic ground states for J 3 < 0. The dark colored spins indicate magnetic unit cell and the numbers are sublattice 
indices. The polar coordinate of the magnetic moments on zth sublattice, uji — {0i,(l)i) is given in Tab. ID 
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TABLE I. Polar coordinate of the moments on different sub¬ 
lattice for the magnetic orders in Eig. [Ha)-(c). 

sublattice configuration with the simple cubic magnetic 
unit cell and zero net magnetization. In Fourier space, 
these states are given by a linear combination of plane 
waves with wave numbers in the set Q = 
with Qx = ( 27 r, 0 , 0 ), q^ = ( 0 , 27 r, 0 ), and q^z = ( 0 , 0 , 27 r) 
(Ref. 1221 ). To investigate the effect of J3, we studied the 
classical ground state by using Luttinger-Tisza method^^ 
In the Luttinger-Tisza method, the ground state is de¬ 
termined by the minimum eigenvalue of 

= ( 4 ) 

i 

where Jij is the 3 x 3 matrix that represents interaction 
between ith and jth spins. The spins are determined by 
a superposition 

Si= ^ (5) 

a=x,y,z 

where a is sum over q^ G Q, and the vector coefficients 
Aa must be further constrained. When J3 = 0 , J(q) is 
diagonal for arbitrary q, with the minimum located at 

q G Q. 

When J3 > 0 , the the coefficients must be chosen so 
that Aa is parallel to q^. The magnetic ground state is 
given by an arbitrary combination of these 3 modes that 
satisfies local constraint jS^j = 1. An example is shown 
in Fig.lHa). Defining polar coordinates 0 and 0 as shown 
in Fig. [I{a), arbitrary 0 G [0,7r] and <p G [—7r,7r] define 
the ground state manifold. It reduces to a two sublattice 
collinear structure for 0 = 0,7r/2 and (j) = n 7 r /2 {n = 
0 , • • • , 3 ), while in general it has four distinct sublattices 


for other sets of {0^(j)). The direction of each spin in 
polar coordinates is given in Tab. [H 

For small J3 < 0 , it is instead required that Aa is 
normal to q^. The magnetic ground state is given by 
an arbitrary combination of these 6 modes that satisfies 
local constraint jS^j = 1 . The solution consists of two 
manifolds of states as shown in Fig. [TJb) and [D^c): a 
double-q state in which two out of the three components 
Aa are finite [Fig. [Hb)], and a triple-q state which all 
three components are finite [Fig. [T][c)]. 

The double-q state consists of two vector components, 
qi and q2, which are elements of Q. Figure [T][b) shows 
an example with qi = q^ and q2 = qx- In this phase, 
the spin axis A2 is fixed parallel/antiparallel to qi, while 
Ai may be in an arbitrary direction. The explicit direc¬ 
tions of spins in the double-q state is given in Table [H 
and an example indicated in Fig. [Hb). In the (6>, 0 ) 
notation, the special values ^ = 0,7r/2 reduce to a single- 
q two-sublattice configuration, while other values with 
0 G (0,7r/2) gives a four-sublattice double q state. 

Figure [T][c) shows an example of the triple-q state. In 
this state, the coefficients A^ are orthogonal to one an¬ 
other. The four sublattice spin pattern collapses to a 
double-q one for 0 = 0,7r/2, and a single-q state for 
(j) = n 7 r /2 {n = 0 , ••• , 3 ). We note that with increas¬ 
ing J3 < 0, the four-sublattice orders eventually become 
unstable. The condition that the above solution is the 
energy minimum gives J2 > l^/sl/S. Hence, if | J3I is suf¬ 
ficiently small, the above-mentioned phases are expected 
to be stable. 

Quantum order by disorder: We next study how the 
quantum fluctuations lift the remaining degeneracy of 
the classical ground state manifold. We calculated the 
quantum correction to the classical ground state energy 
by spin wave analysis. Here, the unit of energy is taken 
as Ji 5 ', where S is the size of spin. The calculation were 
done numerically with Ng = 16 ^ magnetic super cells. 

Figure [ 2 ] shows the quantum correction to the ground 
state energy for J 3 = 1 and J 2 =0.2. 0 and (j) are defined 
as in Fig.[T][a) and TablelH We find that single-q collinear 
antiferromagnetic (AFM) states [6> = 0 or { 0 , (j)) = (f, 0 )] 
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FIG. 2 . Quantum correction to the ground state energy for 
different ground states at Ji = 1 , J 2 = 0 . 2 , J 3 = 1 . 0 and (j) 
are defined as in Fig. [2a). 



FIG. 3. Quantum correction to the ground state energy for 
different ground states calculated by the spin wave analysis at 
Ji = 1, J 2 = 0.2, and J 3 = — 1: (a) double-q and (b) triple-q 
magnetic ground states. 0 and 0 in (a) and (b) are defined as 
in Figs. [TJb) andfljc), respectively. The inset in (a) shows (j) 
dependence of 6E/JiS at ^ = 7 r/ 2 . 

are favored over multi-q states. Specifically, the stag¬ 
gered magnetization in each quantum ground state is ori¬ 
ented along a ( 100 ) axis. This is a natural result as the 
fluctuations tend to favor collinear orders. 

A similar behavior is also found for J3 < 0. Figure 
shows the result for the double-q state [Fig. [ 2 b)]. The 
results indicate that single-q collinear AFM states {0 = 
0,7r/2) are favored over double-q orders. Furthermore, 
fluctuations lift the U(l) degeneracy of the single-q state 
at 0 = 7r/2. The inset in Fig. shows the quan¬ 
tum correction at 0 = it/2] the transverse axis is ^/tt. 
Here, 0 = 0 corresponds to the collinear state with the 


spins pointing along the (100) direction, while they point 
along ( 110 ) for 0 = 7r/4. The minimum is located at 
(^, 0 ) = (7r/2,7r/4), indicating the single-q state with 
q = (0,0, 27 t) is favored by quantum fluctuations; the 
moments point along (110) direction. 

Estimation of the quantum correction for the triple- 
q ground states are shown in Fig. ^h). Similarly, the 
results indicate that, in the triple-q manifold, collinear 
states at 0 = 0 and (6>, 0 ) = (7r/2,7r/4) have the lowest 
ground state energy. This is the ( 100 ) collinear AFM 
state mentioned above. However, all of the states have 
higher energy than the ( 110 ) collinear phase. Hence, 
quantum fluctuations select (110) collinear order. 

The favoring of the ( 110 ) direction for J3 < 0 is a 
highly quantum effect. It is very different from conven¬ 
tional magnetic anisotropy induced by crystal field ef¬ 
fects. In the latter case, a Landau theory analysis in 
powers of the (staggered) magnetization generally ap¬ 
plies, and the leading effect occurs at fourth order. These 
terms favor either (100) or (111), depending on the sign, 
but not (110). 

The above discussion pertains to the direction of the 
staggered magnetization, but we may also ask about 
the uniform susceptibility x(n) (n = h/|/i|) with the 
four-sublattice magnetic structure. Expanding the en¬ 
ergy in small deviations from the ground state, we find 
that x(n) for a given configuration is given by x(n) = 
where Siajp = d'^e/diiiadiiji3 and = 
dm/diiioi with 5 and m being ground state energy and 
magnetization along n, respectively. Here, jj^ia gives the 
small deviation of the ith moment along the orthogonal 
direction a = 1 , 2 . Here, we defined the basis of so 
that £ia,ji 3 become orthogonal, i.e., £ia,ji 3 = £ia^i,j^a,i 3 - 
Using this formula, x(ii) the collinear states are given 
by x(n) = sin^ ^n /(4 + J 3 ), as expected for collinear mag¬ 
netic orders. Here, 0 ^ is the angle between the collinear 
moments and n. This indicates that, for J3 > 0 , the 
maximum of susceptibility will not be along (111). 

To recapitulate the order by disorder analysis, we 
found that, within linear spin wave theory, quantum fluc¬ 
tuations select collinear ground states of the frustrated 
fee antiferromagnet, and in these states the local mo¬ 
ments are oriented along the (100) or (110) axis, depend¬ 
ing upon parameters. The magnetic susceptibility, on the 
other hand, is maximal in {100} planes for J3 > 0. As 
we will see, this behavior contrasts that of the magnetic 
response in the paramagnetic phase, which is favors the 
(111) directions. 

Paramagnetic response: We next study the magnetic 
anisotropy in the paramagnetic phase. Eor this purpose, 
we evaluate the free energy per spin / in an external mag¬ 
netic field. We carry out a high temperature expansion 
of energy / in powers of 07 ^, where 0 is the inverse tem¬ 
perature and i = 1 , 2 , 3 . The expansion of the partition 
function gives 

Z = ^exp - y + 0[(/3Ji)3]|^ (6) 
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where Z/j, = Tr exp(—and {0)h = 

Z^^Tr O exp(—As the contribution to Z from the 
commutator in Eq. [ 6 ]is in the order of 0[{/3Ji)^]^ up to 
the second order in the expansion, we can ignore this 
term. Hence, Z is approximated by 

Z = Zhi^l- p{Ho)h + + 0[{pJif]^ . (7) 

The second order correction to the free energy is given 
by 


/3/(2)(/3h) = (/3Ji)2/f (J2 /Ji,/3/i) 

HPJi? I + gf\lih) y] hi 


(8) 


where, h = |h| is the strength of external field, and ha = 
ha/h is the direction of external field. The contribution 
to the anisotropic term {Ph) of the free energy comes 
only from ([/^Ja (Si ■ %)(Sj ■ (5ij)]2). The explicit 

form of g^\ph) is given by 


9f\m = ^ {S{S + 1) - + 2{S^)^y . (9) 


Here, {S^)h and {S^‘^)h are the thermal averages of 
and 5'^^, respectively, which are functions of the external 
field Ph and the spin size S. Because {^h) > 0 , the 
result indicates that free energy is minimized when the 
external field is applied along ( 111 ). Hence, we obtain 
( 111 ) easy axis magnetic anisotropy regardless of the sign 
of J 3 . We note that this anisotropy is present even for 
S = 1/2. In this case, Eq. ([9|) becomes 


9f\ph) 


= — tanh^(/3h/2). (10) 

5 = 1/2 64 V / y V y 


Hence, for the model in Eq. o, spins favor the (111) 
axis in high temperature. On the other hand, for J 3 > 0 
(J 3 < 0), the ground state is a collinear antiferromagnet 
with spins oriented along the (100) ((110)) axes. The 
distinct difference of the preferred axes at low and high 
temperatures constrasts with the common behavior of 
magnetic anisotropy, which is consistent at all tempera¬ 
tures. This can be understood as a consequence of the 
nature of fluctuations that dominate in different temper¬ 
ature regions. When T ^ Ji, long wave-length fluctua¬ 
tions dominate, lifting the accidental degeneracy of the 
ground states. On the other hand, in the high temper¬ 
ature limit, the spins are strongly disordered, and the 
spin correlation is limited to short range. Hence, in this 
regime, short range correlations between nearby spins is 
most important. Due to the difference in the nature of 
dominant fluctuations, the consequent phenomena can 
be different in the two regimes. 

This understanding also indicates that the information 
on direction-dependent spin interactions is reflected in 
the high temperature paramagnetic phase in an unusual 
way. Hence, this anisotropy-switching phenomena is po¬ 
tentially useful as a probe to experimentally narrow down 
the effective model for heavy element magnets, where var¬ 
ious unconventional interactions can appear due to strong 
spin-orbit interactions. 

Experimentally, a promising candidate to observe 
anisotropy switching is the family of double-perovskite 
compounds with heavy element ions^^“— In these ma¬ 
terials, strong spin-orbit coupling can induce strongly 
anisotropic interactions. Another favorable aspect is the 
suppression of magnetic ordering by geometrical frustra¬ 
tion. As the magnetic anisotropy in the paramagnetic 
phase arises only at higher order in the inverse temper¬ 
ature, a low transition temperature should be favorable 
to enhance these higher order effects. 
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